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Abstract. Introducing sets of constraints, we define new classes of random-matrix 
ensembles, the constrained Gaussian unitary (CGUE) and the deformed Gaussian 
unitary (DGUE) ensembles. The latter interpolate between the CUE and the CGUE. 
We derive a sufficient condition for GUE-type level repulsion to persist in the presence 
of constraints. For special classes of constraints, we extend this approach to the 
orthogonal and to the symplectic ensembles. A generalized Fourier theorem relates 
the spectral properties of the constraining ensembles with those of the constrained 
ones. We find that in the DGUEs, level repulsion always prevails at a sufficiently short 
distance and may be lifted only in the limit of strictly enforced constraints. 
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1. Introduction 

Ever since Wigner introduced random-matrix theory (RMT) in nuclear physics in the 
1950s [1 , that theory has found wide apphcations in modehng the fluctuation properties 
of spectra and wave functions of complex systems ranging from vibrating crystals, to 
microwave resonators, to quantum dots, to atoms, and to the Dirac operator in lattice 
QCD 13 EI- Depending on the symmetry of the Hamiltonian, there are three "classical" 
random-matrix ensembles: the Gaussian orthogonal ensemble (GOE), Gaussian unitary 
ensemble (GUE), and Gaussian symplectic ensemble (GSE) which consist of real 
symmetric, complex Hermitean, and complex quaternion matrices, respectively. The 
corresponding matrix elements are independent, Gaussian-distributed random variables 
with (3 = 1,2, and 4 degrees of freedom, respectively. One of the hallmarks of RMT 
is level repulsion, i.e., the probability of flnding a spacing s between two closely spaced 
neighboring levels is proportional to s^. 

Many complex physical systems exhibit spectral fluctuation properties that are 
in agreement with those of the Gaussian random-matrix ensembles, although the 
corresponding Hamiltonian matrices have structures that differ considerably from a 
Gaussian random matrix. Many-body systems with /c-body interactions, for instance, 
have sparse Hamiltonian matrices, since many-body states that differ in the occupation 
of more than k single-particle orbitals cannot be connected by the interaction. Only 
a few analytical results are known for these fc-body embedded Gaussian ensembles 
|3], and we refer the reader to the recent reviews jSl Ej- Realistic random-matrix 
models for nuclei and atoms also have to include spin and isospin symmetry. The 
resulting two-body random ensemble (TBRE) [71IH] is mathematically very complicated, 
and virtually no analytical results are known regarding the fluctuation properties of 
this important random-matrix model [0]. Another example is given by quasi one- 
dimensional disordered electronic systems. The corresponding Hamiltonians are band 
matrices with zero elements outside a small band around the diagonal. The theoretical 
description of random band matrices is possible due to the simple structure of the 
Hamiltonian |lUj . 

The embedded Gaussian fc-body ensembles, the TBRE, and the ensembles of 
random band matrices have one property in common. They can be viewed as constrained 
random-matrix ensembles. Because of the wide occurrence of such matrix ensembles in 
various branches of theoretical physics, the understanding of their spectral fluctuation 
properties offers a considerable challenge to theory. So far, the available evidence 
(mostly based on numerical simulations for matrices of rather small dimensions) points 
towards fluctuation properties of standard RMT type. The present paper is a flrst step 
toward a common theoretical treatment of constrained ensembles. These ensembles are 
introduced in Sectional and are deflned as Gaussian ensembles of matrices where certain 
linear combinations of matrix elements vanish. We note that the GOE can ultimately 
also be viewed as a constrained ensemble; it is obtained from the GUE through the 
constraint of vanishing imaginary parts of all off-diagonal matrix elements. Remarkably, 
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some properties of constrained Gaussian random-matrix ensembles do not depend on 
the details of the constraints but only on their number and symmetry properties (as 
encoded in systematic degeneracies). In this work, we focus mainly on these rather 
general properties. 

This article is organized as follows. Section El introduces constrained Gaussian 
random matrix ensembles. In Sectional we employ the Harish-Chandra Itzykson Zuber 
(HCIZ) integral formula to express the joint probability function for the eigenvalues 
of the constrained ensemble through an integral over the constraints. In Section |31 
we derive a sufficient condition for the existence of quadratic, GUE-like level repulsion. 
The condition is formulated as a simple inequality that relates the number of constraints 
and the number of degeneracies to the dimension of Hilbert space. We illustrate this 
result with a number of examples. In Section El we analyse the spectal width and 
the distribution of matrix element of the CGUE. In Section IHl we show that the joint 
probability function of the constrained ensemble is related to that of the constraining 
ensemble by a generalized Fourier transform. In Section |3 we consider deformed 
matrix ensembles. These ensembles interpolate between the GUE and the constrained 
ensembles and allow us to study situations where the constraints are slowly switched 
on. We conclude with a Summary. 



2. Constrained Gaussian unitary random— matrix ensembles 

We consider Hermitean matrices acting on a Hilbert space Ti of dimension N. Together 
with Ti, we also consider the linear space V spanned by the Hermitean matrices acting 
on Ti: Every linear combination aA + bB of two such matrices A, B with real coefficients 
a, b is also a Hermitean matrix acting on Ti. In the linear space V, we introduce the 
canonical scalar product in terms of the trace 

{A\B) = Tr(AS) (1) 

for every pair A, B of matrices. This allows us to define an orthonormal basis of iV^ 
Hermitean basis matrices B^ = Bj^ in V which obey 

{B^\Bp) = TiiB^Bp) = (2) 

and 

5: (5,1 = 1, (3) 

where 1 is the unit operator in V. Such a set of basis matrices is, for instance, given 
by matrices that have a unit matrix element somewhere in the main diagonal and 
zeros everywhere else, or an element (1/a/2) somewhere above the main diagonal and 
its mirror image below the main diagonal and zeros everywhere else, or an element 
(i/\/2) somewhere above the main diagonal and its complex conjugate in the mirrored 
position below the main diagonal and zeros everywhere else. This choice of a basis is 
but one example. Any other basis obtained by orthogonal transformations in V (with an 
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orthogonal matrix of dimension N^) from the one just described is equally admissible. 
In general we have no preference in this respect. 

Any Hermitean matrix H acting on Ti. can be expanded in terms of the A^^ 
Hermitean basis matrices as 

H=J2 haB^- (4) 

a=l 

Thus, H can also be viewed as a vector in V. 

We consider the decomposition of V into two orthogonal subspaces labeled V and 
Q. The decomposition is defined in terms of orthogonal projection operators 

Np 

p=l 

Q= E \B,){B,\. (5) 

q=Np+l 

We have 

= V , = Q, V^ = V , Q'' = Q,VQ = 0,V+Q = 1. (6) 

When we sum over matrices in "P-space or in Q-space, we will sometimes use the 
summation indices p and q, respectively, without referring to the partition explicitly 
while a summation over all A^^ basis matrices is indicated by Greek summation indices. 
The operators V and Q are defined in V and have dimension Np and A"q = A^^ — Np, 
respectively. We have been rather explicit in the construction of V and Q because these 
operators differ from the projection operators often used in Hilbert space: The latter 
are defined in terms of the basis vectors spanning Hilbert space, while our projection 
operators are defined in terms of the basis matrices B^ and operate in V. 

In this article, we study constrained ensembles of Gaussian random matrices. Every 
member H of the ensemble is constrained to have zero projection onto Q-space 

QH = or VH = H . (7) 

This condition can also be expressed as 

{By\H) = for allqoT H = J2 Bp{Bp\H) . (8) 

p 

We combine condition (jHJ with the assumption that the ensemble has a Gaussian 
distribution. Using the constraint in the form of the first of Eqs. (jHl), we write the 
probability density W-plH) as 

Wr{H)d[H] = (27r)-^-/2exp {^-l{H\H)^l[S{{B,\H)) d[H] . (9) 

Here d[H] stands for the product of the differentials of the A^^ independent matrix 
elements of H. Eq. © defines the constrained Gaussian ensemble of random matrices 
for the set {Bq} of constraining matrices. 

Every unitary transformation U in Hilbert space induces a transformation Bq 
UBqW of the set of matrices defining Q-space. We say that the set {Bq} is invariant 
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under unitary transformations if the vector spaces spanned by the sets {Bq} and 
{UBgU'^} are identical (This is not the case in general). Let us consider an arbitrary 
unitary transformation H UHU'^ of the matrices in the ensemble (jH)). All terms on 
the right-hand side of Eq. © remain unchanged except for the arguments of the delta 
functions. Here, the unitary transformation is tantamount to replacing Bq UBgW . 
Therefore, W-p is unitarily invariant if and only if the set {Bq} has this invariance 
property. Lack of unitary invariance implies, in general, that the eigenvalues and 
eigenvectors are correlated random variables. 

This is problematic only at first sight. Indeed, the transformation {Bq} — > 
{UBqW} is exactly compensated by the transformation {H = J2phpBp} {UHW = 
J2p hpUBpW}. Instead of considering the sets of constraints {Bq} and {UBqW}, we may 
consider the given set {Bq} and two constrained ensembles, one comprising the matrices 
{H} as defined by Eq. and the other, all matrices obtained from that set by the 
operation H UHU'^, with U fixed. The matrices H and UHU'^ possess identical 
spectra, and in each ensemble the corresponding sets of eigenvalues are encountered 
with equal probability. Therefore, the two ensembles can differ only in the ordering of 
the eigenvalues, i.e., by permutations. Such a difference can arise only if the eigenvalue 
distribution of the original ensemble W-p is not symmetric. 

We group the sets {Bq} of constraining matrices into equivalence classes. Together 
with a given set {Bq}, each class comprises all sets {UBqW} obtained from {Bq} 
by arbitrary unitary transformations U m N dimensions. The "superensemble" W 
is the union of all ensembles W with constraints in the same equivalence class. By 
construction, W is unitarily invariant. We refer to W as to the constrained Gaussian 
unitary random-matrix ensemble (CGUE). The CGUE is ergodic if the members of the 
equivalence class possess symmetric eigenvalue distributions (i.e., distributions that are 
invariant under all permutations of the eigenvalues). In this case, any result derived 
for the superensemble also holds for each of its member ensembles. An example is 
given by the GOE as obtained by constraining the GUE. The GOE is invariant under 
orthogonal transformations, and the eigenvalue distribution must, therefore, be totally 
symmetric. Therefore, the superensemble is ergodic. A counterexample is given by 
the constraint that H have zero non-diagonal matrix elements in the first row and 
column. The eigenvalue distribution factorizes and is not symmetric. The constraint 
is unitarily equivalent, among many others, to the constraint that H have zero non- 
diagonal matrix elements in the fcth row and column, with k = 2, . . . , N. Thus, results 
derived for the spectral fluctuation properties of the associated CGUE apply jointly 
to all these constrained ensembles but not to the individual members. We see that 
our focus on the CGUE prohibits the study of unsymmetrical eigenvalue distribution 
functions which might be useful, e.g., in studies of symmetry-breaking. It may happen 
that the vector space spanned by the set {Bq} is invariant under a subgroup of the 
unitary group. In that case, the equivalence class is generated by the coset space of the 
unitary group with respect to that subgroup. 
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For a set of constraints {Bg}, the CGUE is defined by 
Wr{H)d[H] = 

(27r)-^-/2exp [-^-{H\H)y[H] J d[f/](n 6{{UBgU^H))^ . (10) 

The integral d[U] extends over the unitary group in N dimensions. The ensemble is 
obviously invariant under unitary transformations of H. The Haar measure of the 
unitary group is normalized to one, i.e. 



dp] = 1 . (11) 

Proceeding as usual, we diagonalize the matrix H with the help of a unitary matrix V, 

H = VxV^ , (12) 

where x = diag(xi, . . . , xn) is the diagonal matrix of the eigenvalues. Here, and in what 
follows, we denote diagonal matrices by small letters. The integration measure becomes 

d[H] = ^ A\x)d[x]d[V] , (13) 

where A(x) denotes the Vandermonde determinant 

A(x) = Yl (xk- X,) . (14) 

Eq. (jl3|) shows that eigenvalues and eigenvectors are uncorrelated random variables. 
The distribution of the eigenvectors is defined by the Haar measure. Therefore, the 
joint probability distribution Pp{x) of the eigenvalues is the object of central interest in 
this paper. It is given by 

(27r)(^«-^)/2 



Pvix) = M^l^exp (-i(x|x))A2(x)Fp(x) , (15) 



where 



FAx)^ I d[U][^l[5i{Bg\UxU^))'^ . (16) 



This construction of the CGUE may seem unnecessarily involved. Why not start 
from the constrained ensemble defined in Eq. use the transformation of variables 
and measures as given by Eq. (|T3|) . and integrate W over the unitary group to obtain 
the distribution function for the eigenvalues? Following this path, we actually arrive at 
Eqs. flT3j) and (fT^ and, thus, at the CGUE: In this approach, the eigenvalue distributions 
are always symmetric under permutations of the eigenvalues, since the integration over 
the unitary group sums over all N\ different unitary matrices that diagonalize a given 
matrix H with non-degenerate eigenvalues. We believe that our line of reasoning 
shows more clearly the conceptual framework we use. 

It would be desirable to have a classification scheme for the CGUEs. How many 
such ensembles are there for a given number of constraints? Equivalently, how many 
different equivalence classes exist for a fixed number of constraining matrices? We have 
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not investigated these questions yet. In the present paper, we focus attention on joint 
properties of the CGUEs. 

Because of the presence of the function F-p[x) on the right-hand side of Eq. (fTHjl . 
the spectral statistics of the CGUE differ from those of the GUE. To work out that 
difference, we must study the function -Fp(x). 



3. The HCIZ Integral 

To calculate F-p^x), we replace the delta functions by Fourier integrals, introducing the 
set {tq} of Nq integration variables. Then, 

Fp{x) = (27r)-^« Jli^tgJ d[f/]exp (^zJ2t,{Bq\UxU^)^ . (17) 

The integral over the unitary group is the Harish-Chandra Itzykson Zuber (HCIZ) 
integral [TTl [Tallin]- It can be expressed in terms of the eigenvalues bj{{tg}), j = 1, . . . ,N 
of the matrix 

Bi{t,}) = J2t,B,, (18) 
and is given by 

The determinant has Xk in the fcth row, and 6/ in the Ith column, with k,l = 1, . . . ,N. 
We thus obtain 

p _ nfji' j! fyr det [exp jtXkbi)] 

The function F-p{x) is obviously symmetric with respect to all permutations of the 
eigenvalues Xk- The integrand of F-p{x) in Eq. piU) has that same property with respect 
to the eigenvalues bi of the matrix B. Both properties are a consequence of the unitary 
invariance of the CGUE. Multiplying every eigenvalue Xk by a parameter r so that 
Xk — > TXk for all k, we get 

Fv{tx) = r^^«F-p(x) . (21) 

This shows that F'p{x) is a homogeneous function of the eigenvalues which has a 
singularity of order Nq when all eigenvalues tend to zero simultaneously. However, 
F-p (x) does not diverge generically when any one of the eigenvalues vanishes individually. 

Using Eq. in Eq. (fTSj) . we find for the joint probability distribution of the 
eigenvalues the expression 

p _ exp(-i(x|3:))A^(x) r det[exp (z&.x^)] 



^2t,){N+Nq)/2 ^N{N-l)/2^\ J 11 A(6)A(x 
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4. Sufficient condition for level repulsion 



4.I. Unitary case 

We use the function F-p{x), as given in Eq. (plj) . to investigate the spectral statistics 
and, in particular, level repulsion at short distances for the constrained ensembles. It 
is obvious that GUE-like level repulsion, as given by the factor A^(a;) in Eq. 
will prevail unless the function F-p{x) has singularities whenever two eigenvalues Xj,Xk 
coincide. This prompts us to study especially the singularities of F-p{x). 

The integrand in Eq. (j^UI) has no singularities: Any zeros of the denominator are 
canceled by corresponding zeros of the numerator. Indeed, zeros of the denominator 
arise when two or more eigenvalues bj of the matrix B coincide. We assume that 
61 = 62 = • • • = are L degenerate eigenvalues. Without loss of generality, we may 
assume that they vanish. (Under translation bj bj + Ab with fixed Ab for all j, 
the integrand is multiplied by the non-singular function exp {iAbJ^^k))- Then the 
integrand becomes 



det[exp (ibjXk)] 



A{x)A{b) 



det 



6l=0 

1 
1 



(23) 



n 



Xi 
X2 



X2 



jixibL+i 
1 pix2bL+i 



1 Xn 



X^ 



^ix\bj^ 



^ix]\[b]\[ 



l\ 



lln=LH 



A(5l+i, . . . ,6Ar)A(x) 



(where A(6i+i, . . . , ^at) denotes the Vandermonde determinant of the eigenvalues 
6^+1, . . . ,67V) and is obviously finite. The same conclusion applies when two or more 
eigenvalues Xk coincide. For a discussion and alternative calculation of the HCIZ integral 
in the presence of degeneracies, we refer the reader to Refs. jl4[ I15j. 

We conclude that the A'^g-dimensional integral in Eq. (j^Ull can diverge only because 
the domain of integration is not bounded. This suggests that we introduce Nq- 
dimensional spherical coordinates Ilgdtg = t'^'^'^dtdVL (with t the radial variable and Vt 
the set of angular variables), and focus attention on the radial integral 



dt 



_j^det[exp {ibjXk)] 



(24) 



From Eq. (jl8|) . it follows that the eigenvalues are linear homogeneous functions of the 
variables tq so that bj{{Xtq\) = \bj{{tg}). This implies bj(t,Q) = tbj{l,Q). 

We distinguish two cases, (i) There are no systematic degeneracies among the 
eigenvalues bj, i.e., there exists no domain of integration with measure > in which 
eigenvalues are degenerate. (Accidental degeneracies that occur on a set of measure 
zero of the integration domain are irrelevant.) Then, the Vandermonde determinant in 
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the denominator yields a factor and the radial integral becomes 

oo 

J dt t^'3-^-^(^-^)/Met[exp(zt6j(l,r])a;fc)] . 



(25) 



The generalization of Eq. ()23|1 to L = N degenerate eigenvalues shows that the integrand 
is not singular at t = 0. The radial integral is sure to converge if Nq < N{N — l)/2. (ii) 
There is a set of L systematically degenerate eigenvalues. The domain of degeneracy 
extends out tot = oo. In this case, we use the integrand in the form ()2Hj] . For the radial 
integral, this yields 



^NQ-l+L{L-l)/2-NiN-l)/2 



Xi 
X2 



X-^ 
X2 



„L-l 



^ix2bN 



^ixj\fbj\f 



(26) 



i Xjv • • • 

Again, the integrand is well behaved at t = 0, and the integral is guaranteed to converge 
for Nq + L{L — l)/2 < A^(A^ — l)/2. This result can be generalized. Let there be J sets 
of Lj,j = 1, . . . , J degenerate eigenvalues with degeneracy domains that extend out to 
t = 00. Then, quadratic GUE-like level repulsion is guaranteed to occur if the number 
of constraints Nq, the numbers of degenerate levels Lj, and the matrix dimension N 
obey the inequality 

J 



iVQ + ^L,(L,-l)/2<iV(iV-l)/2. 



(27) 



The condition (j2Zj) is a sufficient (but not neccessary) condition for GUE-type level 
repulsion to occur in the presence of Nq constraints. (It may happen that the 
exponential functions in the integrand provide convergence even if condition ()27|1 is 
violated.) The inequality ^T7\i makes no reference to the specific structure of the 
constraints; the only input is the number of constraints and the number of degeneracies. 
We recall that degeneracies are related to symmetries. In this sense, the number and 
symmetry properties of the constraints are at the root of the inequality ()27p. 

The possibility to derive the general result (fTfjl without being specific about the 
nature of the constraints is, of course, due to the unitary symmetry of the CGUE. 
We have mentioned already that the function F-p{x) will affect spectral fiuctuation 
properties not only on a very short scale (level repulsion) but also on larger scales. It 
appears, however, that that infiuence depends on specific properties of the constraints 
and requires explicit calculation of F'p{x). 



4.2. Orthogonal and symplectic case 

Instead of considering the Hilbert space of Hermitean matrices, we might have started 
out from the Hilbert spaces of real symmetric or complex-quaternion matrices. In both 
cases, we follow the same steps as for the CGUE and arrive at the analogue of Eq. (fTH|) 

P-p{x) oc exp (^~^{x\x)^\A{x)fF-p{x) . (28) 
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The factor |A(a;)|'^ stems from the Jacobian of the eigenvalue/eigenvector transformation 
for real symmetric (/5 = 1), and complex-quaternion matrices {jS = 4). We have omitted 
the overall normalization constant. Likewise, we have not distinguished the functions 
P-p and F-p from their unitary counterparts as there is no room for confusion. The 
function F-p is analogous to Eq. (fTTjl and reads 

Fp{x) oc J Y[dtgJ d[Uf,]exp (^iY,t,{Bg\UfsxUl)^ . (29) 

We denote orthogonal and symplectic matrices as Ui and f/4, respectively. The 
integral dlUp] is over the orthogonal (/5 = 1) or symplectic (/3 = 4) group in N 
dimensions. The constraining matrices Bg are real and symmetric for (3 = 1 and 
complex quaternion for /3 = 4. Eqs. (|^ and define the spectral statistics for 
the constrained Gaussian orthogonal ensemble (CGOE) and the constrained Gaussian 
symplectic ensemble (CGSE). 

Progress is hampered by the fact that the analogue of the HCIZ integral is not 
available in closed form for the orthogonal and the symplectic groups. However, we 
might employ an asymptotic expansion of the HCIZ integral that has recently been 
re-derived in Ref. ^H]- We assume that neither the eigenvalues x nor the eigenvalues 
of the matrix B possess any systematic degeneracies. Again, we introduce spherical 
coordinates Ug dtg = t^Q'^dtdn. We are interested in the behavior of the integrand 
for large values of the radial coordinate t. Since bi(t,Q) = tbi{l,Q), the differences of 
eigenvalues of B grow linearly with t, and we may confine ourselves to the leading term 
in the asymptotic expansion 

/ d|C/,l exp (ij:t,{BM-Ul)) « i^^lTAmi?/^' • 



q 



Correction terms are inversely proportional to powers of the differences of the eigenvalues 
{bj — bk). Following the same arguments as in the unitary case, we find that the radial 
t-integration is guaranteed to converge if 

NQ<f3N{N-l)/A. (31) 



In the absence of any systematic degeneracies of the constraints, the inequality ()3Tjl 
is a sufficient condition for the existence of eigenvalue repulsion of the canonical form 
\xi — Xj\^ in the constrained Gaussian ensembles with Dyson index (3. The more general 
result (|?fjl for /3 = 2 and the inequality (pT|l are central results of our work. 



4-3. Examples 

We present several examples for the CGUE, the CGSE, and the CGOE, which illustrate 
the power and limitations of the inequalities derived above. 

(i) We consider the GOE as obtained by constraining the CUE. The P-space 
consists of real symmetric matrices and the Q-space of purely imaginary antisymmetric 
matrices. We have Nq = N{N — l)/2 and no systematic degeneracies. The 
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condition (jTTj) is violated, as must be the case, because the eigenvalue distribution 
function of the GOE carries |A(x)| in the first power. This implies that 

oc ^ . (32) 

This relation can be proved by employing a relation between complementary distribution 
functions that is derived in Section IHl see also the Appendix. 

(ii) We consider the case where Q-space is spanned by the set of diagonal matrices. 
The number of constraints is and the condition ()27|) is fulfilled. Level repulsion 
prevails. For A^ ^ 1, the corrections to the GUE average level density and the two-point 
function can be worked out analytically. We use the supersymmetry technique and the 
standard saddle-point approximation. This approximation consists in neglecting terms 
which are small of order 1/A^ and results in an asymptotic expansion in inverse powers 
of A^. The deviations from the GUE caused by the constraints are of order 1/A^ and, 
therefore, do not affect the saddle-point condition. Hence, we find that the corrections 
to the GUE results vanish like inverse powers of A^. 

(iii) This case is converse to the previous one: P-space is spanned by the set of 
diagonal matrices. The number of constraints is Nq = — N, and condition ()27j) 
is violated. The levels in P-space obey Poisson statistics. Therefore, the function F-p 
must compensate the factor A^(x) in Eq. (j221)- Hence we must have 

^ • (33) 

It seems difficult to prove the relation ()33|) by direct calculation. 

(iv) Block-diagonal Hermitean matrices: We require that two rectangular blocks 
in H carry zero matrix elements. If the two blocks consist of all non-diagonal matrix 
elements in row k and in column k, the number of constraints is 2(A^ — 1). In this 
case, A^ — 2 of the eigenvalues of B are degenerate (i.e., vanish), and the remaining two 
eigenvalues differ only in sign. The condition (P7j) is violated (albeit weakly). These 
considerations are easily extended to bigger blocks. 

(v) Another example is furnished by the EGUE, the embedded Gaussian unitary 
ensemble with fc-body interactions [H IHl IE] • This ensemble differs from a CGUE in that 
its matrices Bp do not obey the normalization (j2)). However, the corresponding CGUE 
is dense in the EGUE. With / the number of degenerate single-particle states, and m 
the number of identical Fermions, we have A^ = The number of independent k- 

body matrices Bp is Np = 2(Q,^)^ — (J^. According to condition (j7fj) . level repulsion is 
guaranteed if there are no degeneracies and if A^(A^ + l)/2 < Np, or if 

-l) ■ (34) 

For arbitrary /, this condition is met only for k = m, and there it is expected. 
Unfortunately, we cannot draw any non-trivial conclusions for the EGUE. 

(vi) "P-space consists of all traceless Hermitean matrices. The single constraint 
(A"q = 1) is proportional to the unit matrix, and all A^ eigenvalues of the constraint are 
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degenerate. The probability distribution can be worked out by means of Eq. (j2Sl), 
and the result differs from the GUE eigenvalue distribution only by an overall factor 
6{xi + xi + ■ ■ ■ + xn), which yields a GUE spectrum with vanishing centroid. The 
spectral fluctuations are identical to the GUE though inequality ()27|) is not fulfilled. This 
demonstrates that the conditions (P7j) and, by the same token, (p?T|l are only sufficient: 
Their fulfillment guarantees level repulsion with power (3, and they cannot be fulfilled 
if there is no level repulsion with power (3. However, no conclusion can be drawn about 
level repulsion if these inequalities are violated. 

(vii) We consider the GUE as obtained by constraining the GSE. The number of 
constraints is Nq = N{N — 1), there are no systematic degeneracies, and the Dyson 
index is /3 = 4. The inequality is violated, as it must be. 

(viii) We consider an ensemble of real symmetric block-diagonal matrices consisting 
of two blocks with equal dimensions. Clearly, eigenvalues belonging to different blocks 
are uncorrelated. The Q-space is given by the chiral GOE with quadratic blocks. The 
number of constraints imposed on the GOE is Nq = there are no systematic 
degeneracies, and (3 = 1. The inequality is violated, as must be the case. 

(ix) We consider the chiral GOE as obtained by constraining the GOE. Matrices 
of the chiral GOE consist of two rectangular off-diagonal blocks of size n x m. For 
m > n, there are m — n zero eigenvalues, and the remaining 2n eigenvalues come in n 
pairs with opposite signs. The correlation of the positive eigenvalues is determined by 
the expression [TH] 

n 

n i^^-^?in^i'"""'- (35) 

We consider two cases. First, we set m = n. This case is converse to example (viii), 
as the roles of "P-space and Q-space are interchanged. According to expression (jH^j) . 
eigenvalues with equal signs repel each other linearly, but there is no repulsion between 
the smallest positive eigenvalue and the largest negative eigenvalue. The number of 
constraints is Nq = n{n + 1) = {N/2){N/2 + 1), there are no systematic degeneracies, 
and (3 = 1. The inequality (jHT|) is violated, as it must be. For the second case of the chiral 
GOE, we set m = n + 1. In this case, the eigenvalue spectrum has one zero eigenvalue, 
and n pairs of eigenvalues with opposite signs. We have Nq = (n + 1)^ = (A^ + 1)^/4, 
there are no systematic degeneracies, and P = 1. The inequality is violated. 
However, inspection of expression shows that there is linear level repulsion between 
any pair of neighboring levels. In particular, the zero eigenvalue repels the smallest 
positive eigenvalue. Again, we cannot draw any conclusion about level repulsion from 
the inequality (jSH) since it is violated. 

5. Further properties of the CGUE 

Constraints affect not only the spectral fluctuations of the CGUE but also the spectral 
width, and the distribution of matrix elements of the CGUE. To show this for the 
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spectral width, we recall the usual normalization condition H'^^ = X'^/N of the matrix 
elements if^^^ of the GUE. Here N is the matrix dimension, the overbar denotes the 
ensemble average, and 2A is the radius of Wigner's semicircle. With this normalzation, 
the spectral width of the GUE has the value ( Trace if ^ = A^. To work out 
(1/iV) Trace if ^ for the CGUE, we start from the normalized distribution function of 
the eigenvalues, 

Pcgue(x) = ^^A'{x)F{x) eM-^E^') ■ (36) 

Here A/" is a normalization factor. For F{x) = 1 (no constraints) this function yields the 
GUE with the above-mentioned normalization. To work out the normalization factor 
Af, we relate PcGUE(a^) to the distribution function P-p(x) defined in Eq. (fT^ with 
normalization factor A/q = (27r)'^^'5~^^/^/ H k\. We substitute xj = i/jX/ \/N and use the 
fact that F{x) = {y/N/\)^Q to obtain 

Ar=(J)(^^-^^)/^Aro. (37) 

A straightforward calculation then yields 

(lTrace:^)c^UE = A2[l-^]. (38) 

This shows that in comparison with the GUE, the spectral width of the CGUE is reduced 
by the factor \Jl — Nq/N'^. Every constraint - irrespective of whether it is located on 
the main diagonal or not, or whether it affects the the real or the imaginary elements 
of H - yields the same contribution to the argument of the square root. 

For the GUE, the matrix elements are Gaussian-distributed random variables. 
Because of the constraints, this is not so for the CGUE. To show this, we calculate 
the variance 

~ j2 r 1 ^2 



Vogue 



— Trace H'^ 
N 



— Trace if^ 



(39) 



of (l/A^)Trace iJ^ in two ways. Using Pcgue{x) as given by Eq. (j^HI) "we obtain 



W = 1^(1 - ^) . (40) 



On the other hand, Eq. (j38|l implies for the CGUE 



H,.H,, = -(1 - ^) . (41) 



Assuming now that the matrix elements do have a Gaussian distribution, we find for 
the variance of (1/A^) Trace the value 



• ^ ^ 

This disagrees with Eq. pUj) and shows that the deviations from the Gaussian 
distribution increase monotonically with increasing number of constraints. While the 
matrix elements of the constrained ensembles defined by Eq. do have a Gaussian 
distribution, this is not the case for their unitarily invariant counterparts defined in 
Eqs. fjT^ and (fTHj). 
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6. Complementary distribution functions 

Equation ()22|1 expresses the joint probability function P-p{x) of the CGUE as an integral 
over the constraints {Bg}. We can interchange the roles of the P-space and the Q-space 
and study the distribution function Pq using the set {Bp} in P-space as constraints. 
The probability distribution Pq{x) has the form of Eq. with 

Fq{x) = J d[U]^Y[6{{Bp\UxU^))^ (43) 

playing the role of F-p{x). The integrands of either function depend only on the 
eigenvalues b of the matrix B of the constraining ensemble. This suggests that it should 
be possible to express Pp{x) as an integral over the complementary distribution function 
Pq(6), and vice versa. In this Section we show that this is indeed the case and can be 
done via a generalization of Fourier's theorem. 
We start from the identity 

1[6{{B,\H)) = JU^tiJ lidtpexp (^z{J2tMH?j l[^{tp) . (44) 

The matrices B„ form a basis of V. Therefore, the matrix B = J2a taPa is the general 
Hermitean matrix in dimensions. By a suitable orthogonal transformation, we can 
replace the (general) basis matrices B^j by the special basis set defined below Eq. (jH}. 
For the transformed integration variables i^-, we have Ilcrdtcr = rio-dto-. The product 
Her dtff is the product of the differentials of all matrix elements and, hence, equal to 
d[B]. The delta functions S{tp) can be written identically as S{{Bp\B)). Thus, 

l[Si{B,\H)) = f d[B] exp{^{B\H))l[6{{Bp\B)) . (45) 

The identity ()45|) relates the constraints in P-space and in Q-space via Fourier 
transformation. The Fourier integral is taken over Hermitean matrices in dimensions. 

We recall H = UxU\ write B as B = VbV'^ , and integrate the identity over 
the unitary group U, using the definition |TBj) of -Fp(x) and the HCIZ integral. This 
yields 

Fv{x) = J d[U] J d[B] exp {z{B\UxU^))l[6{{Bp\B)) 



(27r)^(^'i)/2-^Q 



^\{5{{Bp\VbV^)) 



d[U] / d[V] / d[b] A\b) exp U{VbV^UxU^) 



^2^^NiN-m-N, det[exp{tb,x,)] 

Nl ,7V(7V-l)/2 J ^i'^J "^y^) A(x) 



X fd[V] USiiBAVbV^)) . (46) 
•' p 

Combining this result with the definition ()43|) of Fq{x), we obtain the identity 

Hx)Mx) = ^ ^; J d[b] det[exp (tb.xi)] A(6)Fq(6) . (47) 
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This identity shows that A{x)F'p{x) is the Fourier transform of A(6)Fq(6). It is a 
generahzation of Fourier's theorem to a class of symmetric functions of the variables 
xi,X2, . . . ,xn- It is the main result of this Section. We recall that F-p{x) and -Fq(6) 
are directly related to the joint probability functions Pr{x) and PQ^b), respectively, via 
Eq. Thus, the complementary probability functions themselves are also obtained 

from each other by a generalized Fourier transform. 

Writing the identity (jTrf) for both F-p{X) and FQ^b), using the reality of F-p and 
Fq in the latter relation, and inserting the result into the former, we find 



A(x)Fp(x) = ^— ^ / d[b] / d[y]det[exp{zxjbi)]det[exp{-ib0i)]A{y)Fr{y) . 



holds for a class of symmetric functions. 

To illuminate the result (@7j), we discuss two examples. In the first example, the 
P-space consists of all real symmetric matrices, and the Q-space of all antisymmetric 
matrices with purely imaginary matrix elements. The joint distribution function 
of the GOE is known, and so is that of the (purely imaginary) antisymmetric 
Gaussian ensemble Both these functions are related via the generalized Fourier 

transform (jTFj) . and the explicit calculation is presented in the Appendix. In the 
second example, the P-space consists of two block-diagonal matrices with equal 
dimensions. Thus, the eigenvalue distribution is that of a superposition of two GUEs. 
The complementary Q-space is given by the chiral GUE used to describe universal 
phenomena of the QCD Dirac operator. The eigenvalue distribution for this ensemble 
is also known [18 . Both eigenvalue distributions are related to each other by the 
generalized Fourier transform. 

7. Deformed Gaussian ensembles 

So far, we have formulated the constraints in terms of delta functions (see Eq. Q). 
It is tempting to consider a more general case where the matrices of the V- and the 
Q-space both have non- vanishing probability but differ in the widths of their Gaussian 
distribution functions. By changing the width of the distribution in Q-space, we are 
able, for instance, to study the transition from the GUE to the CGUE. We refer to such 
ensembles as deformed Gaussian ensembles. 

With {Ba} an orthonormal basis in V, we consider an ensemble of Hermitean 
matrices 




(48) 



Eq. (gHl) implies that 
(iV!)2 



d[b] det[exp (ixjbi)] det[exp {-ibjyi)] = ]J (5(x; - yi) (49) 



1=1 



1 ' 



(50) 
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where Sp and Sg are real independent Gaussian variables with zero mean and variances 

4 = ^^^, = ^- (51) 

Here and in what follows, we indicate ensemble averages by a bar over the quantity 
of interest. We note that for Xp = 1 and Aq = oo, we recover the constrained 
ensemble ©, while for A-p = 1 = Aq, we recover the GUE. We note that the deformed 
Gaussian ensemble fl5()|l is not unitarily invariant unless A^^ = Aq. In calculating the 
spectral distribution function, we do not proceed as in the main part of Section 
where we constructed first the unitarily invariant ensemble. Rather, we use the shortcut 
described below Eq. (fTH|) and integrate directly over the unitary group. This procedure 
automatically generates the deformed Gaussian unitary ensemble (DGUE). 

We make use of = {H\B„), d[H] = Ho- ds^-, and start from the probability density 
formula 

WpQ{H,Xv,XQ)d[H] = (52) 



XA ' Xq\ ' ( Ap^,.._,, Ac 
2ti] \2tx , 



Nr. 

We employ the completeness and orthonormality of {-Bo-}, substitute 

XvY.{H\Bpf + XQY,{H\B,f = Xv{H\H) + (Aq - Ap) ^(^1^-?)' > (53) 

PI Q 

express H in the diagonalized form, H = UxU\ and write d[H] as in Eq. (|T!?|) . Inspecting 
the resulting formula, we find that the distribution function P-pq(x, A-p, Aq) of the 
eigenvalues of the DGUE is given by the integral over the unitary group in dimensions 



^pq(x,Ap,Aq)= — -f 'jv exp(- — (a:|a:))A^(x) 



X 



Np Nq 

-/Aq\^ 

m exp Y.{UxU^W] ■ (54) 



\ I? 
We use the Hubbard-Stratonovich transformation 



exp 



Aq — X-p 



iJJxU^Bqf 



^ J dtqexp (^ttg{UxU^Bg) 



X 
2 



K I (55) 



where 



A = (56) 

Aq — A-p 

and integrate over the unitary group with the help of the HCIZ formula p9|) . We obtain 

o . ^ ^ ^ / Ap\ ^ / Aq\ ^ (27r)^(^-i)/2 / Ap , , .\ , 
PpQix, Xp, Aq) =i^-j ^—j ^, ^^(^_,)/, exp ^--{x\x) j A (x) 

/ A \ 2 /• / X 2] det [exp (ixkbi)] 
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The {hi} are the eigenvalues of the matrix B = Y,q iqBq. We note that for \q ^ X-p = 1, 
we have A ^ oo, the Gaussian factors turn into deha functions, the integral over the 
{tq} becomes a constant, and Ppq{x, X-p, Xq) turns into the GUE distribution. On the 
other hand, for X-p = 1 and Xq — ^ oo, we have A — >• 0, the Gaussian cutoff becomes 
irrelevant, and the right-hand side of Eq. (jFTj) approaches that of Eq. ()22j) . Thus, the 
DGUE correctly interpolates between the GUE and the CGUE. 

We recall the discussion in Section |3] and observe that the integrals on the right- 
hand side of Eq. (|K7jl converge for all positive values of A. This shows that in the DGUE, 
Ppq{x, Xp, Xq) vanishes whenever two eigenvalues Xk,xi coincide no matter how large 
Xq, and level repulsion prevails. If level repulsion disappears, it does so abruptly at 
Xq = oo. On physical grounds we expect, of course, that as Xq increases, level repulsion 
affects ever smaller distances between levels. 

The last term on the right-hand side of Eq. ()57j) can be reinterpreted as the 
ensemble average of det[exp {iXjbk)]/{A{x)A{b)) over the Q-space ensemble of matrices 
B = J2qtqBq with Gaussiau-distrlbuted tq with the variances = 1/A. We have 

Tldt Pxnl^ A 2\ det [exp jtxkbi)] _ r det [exp (zo: A)] 

11 dt, exp ^--t^j ^(^)^(^) - J d[b\PQ{b,X) ^^^^^^^^ (58) 

where (cf. Eq. (fT3jl ) 

(27r)^(^-i)/2 / A \ ^ fx 




PQ{b,X) = ^ i^-j exp [--{b\b)jA\b) J d[U]llS{{B,\UbU^)) (59) 

denotes the disribution function of eigenvalues bk of this ensemble. This shows that 
Ppq{x, Xp, Xq) can be expressed entirely through the rescaled eigenvalue distribution 
PaibA)- 

o . ^ ^ ^ / Ap\ * / Aq\ ^ (27r)^(^-i)/2 / Xp . , .\ . 
Ppq{x,Xp,Xq)= ^—j ^—j -j;^^-^^^^exp\^- — {x\x)jA (x) 

./d„P,,M)^^^^. (60) 

Throughout the derivations in this Section we may, of course, exchange the roles of 
the "P-space and of the Q-space. The resulting equations are obtained by substituting 
P ^ Q except that the term in the exponent on the left-hand side of Eq. is 
now positive. This causes the imaginary unit i to disappear in that and all following 
equations. We note that the generalized Fourier transformation of Section IHl can be 
extended to the DGUE. We also note that the transition from GUE to GOE, or from 
GUE to the antisymmetric Gaussian ensemble, can be described in the framework of 
the DGUE. In this case, one recovers the result given in Ref. |I9j . 

8. Summary 



Introducing sets of constraints, we have defined new classes of random-matrix ensembles, 
the constrained Gaussian unitary (CGUE), and the deformed Gaussian unitary (DGUE) 
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ensembles. The CGUEs consist of Hermitean random matrices where the constraints 
set certain hnear combinations of matrix elements to zero. The DGUEs interpolate 
between the GUE and the corresponding CGUE. Using the Harish-Chandra Itzykson 
Zuber integral for the CGUE, we have found a sufficient condition for GUE-type level 
repulsion to persist in the presence of constraints. This condition depends only on the 
number of constraints and on the internal symmetries of the constraining matrices and 
can be formulated as a simple inequality. We have derived a generalized Fourier theorem 
which relates the spectral properties of the constraining ensembles with those of the 
constrained ones. We have shown that in the DGUEs, level repulsion always prevails 
at sufficiently short distances. It is only in the limit of strictly enforced constraints 
that level repulsion may be lifted. As a result, we find that GUE-type level repulsion is 
remarkably robust. The extension of this approach to the orthogonal and the symplectic 
cases is hampered by the fact that the analogues of the Harish-Chandra Itzykson Zuber 
integral formula are not known in closed form. In these cases, we found a sufficient 
condition for the canonical level repulsion only for certain classes of constraints. 
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Appendix: Relation between the GOE and the antisymmetric Gaussian 
ensemble 

The Gaussian antisymmetric ensemble and the GOE are complementary ensembles. In 
this appendix, we show that the eigenvalue distributions of GOE and the antisymmetric 
Gaussian ensemble are related to each other by the generalized Fourier transform (07j). 
For simplicity, we consider matrix ensembles with matrices of even dimension = 2m. 
We start by considering the GOE basis matrices as the constraints and use the 
probability distribution of the antisymmetric Gaussian ensemble as given in Ref. J7] 

PAb)^Sexp(--m) n ib'.-birUm + bm^i)- (A.l) 

^ ^ l<j<k<m 1=1 

We have omitted an overall normalization constant. The symbol S denotes 
symmetrization with respect to all eigenvalues. According to Eq. one thus finds 

m 

mMb)^A n {b]-hifx{5{h+h^^i)/m (a.2) 

l<j<k<m 1=1 

o^^ll 1 • 

k=l 
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The symbol A denotes antisymmetrization with respect to all eigenvalues. We have 
used the identity 

n i4-xirf[6{x,+x,^J = A{x)f[^-^^^±^, (A.3) 

l<fc</<m j=l j=l -^j -^j+m 

which can proved by expanding the products and exploiting the 5-functions. We insert 
the last relation ()A.2|) into the Fourier relation (P7j) and find 



TT -^det [exp {ihjXk)] 
1=1 



(A.4) 



We have performed one-half of the integrations by means of the (5-functions, and 
bm+j = —bj thus holds in the exponential. We have dropped the antisymmetrizer A 
since the determinant is antisymmetric anyway. We use the Laplace expansion of the 
determinant and obtain terms of the form 
/■ dfe 

n / IT exp (ibkixj^ - Xj^^j) = (ivr)'" J] sign(xj, - xj^^J . (A.5) 
k=i-' fc=i 

Summing up all terms, one finds 

m 

A{x)Fq{x) (xAl[ sign(a;j-, - xj^^J . (A.6) 
fe=i 

This expression is totally antisymmetric under permutations of any two eigenvalues and 
is a homogeneous function of degree zero. Thus, 

A{x)Fq{x) oc sign(A(x)) , (A.7) 

and this yields the GOE distribution once Eq. (fT3jl is employed. 

We can also employ the Fourier transform to go from the GOE to the antisymmetric 
Gaussian ensemble. For this purpose, we insert the GOE expression ()A.7|) into the 
Fourier relation (j47|) and find 

A{x)F-p{x) (X J d[b]sign{A{x))det [exp (ibjXk)] . (A. 8) 

This integral can be performed by using, for instance, the results presented in Sect. 3 
of Ref. HH]. We find 

A{x)Fr{x) oc PfKfc]j,fc=i...2m ■ (A.9) 

Here, 



ajk= / dwdt; (e*^^ V^'=^ - e'^'^V^^'' (A.IO) 



= 8ni{xk — Xj) ^6{xj + Xk) , 
and Pf denotes the Pfaffian 

Pf K'fc] = E(-l)^%ii2ai3i4 • • • aj,^,j2m , (A.ll) 
p 

with P running over all (2m)!/(2™'m!) permutations (ji, . . . ,j2m) of (1, 2, ... , 2m) with 
restrictions ji < ja, js < ^4, • • • j2mi < j2m and ii < < . . . 22m-i- 
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The expression fjA.9|) is thus completely antisymmetric under permutations of any 
two of its arguments. The probability density we seek is essentially this expression, 
multiplied by the Vandermonde determinant A(x), and is thus completely symmetric 
under permutations of its arguments. We may therefore focus on just one term of the 
Pfaffian and write 

Pvix) OC S exp i--{x\x)]A{x) Yl{Xk - Xk+rn)~^5{Xk + Xk+m) ■ (A.12) 

^ ^ fc=l 

Employing the 5-functions and the identity ()A.3|) . we recover Eq. ()A.1|) . 
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